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We consider the interaction of a ferromagnetic spinor Bose-Einstein condensate with a magnetic 
held gradient. The magnetic held gradient realizes a spin-position coupling that explicitly breaks 
time-reversal symmetry T and space parity V, but preserves the combined VT symmetry. We 
observe using numerical simulations, a phase transition spontaneously breaking this remaining sym¬ 
metry. The transition to a low-gradient phase, in which gradient effects are frozen out by the 
ferromagnetic interaction, suggests the possibility of high-coherence magnetic sensors unaffected by 
gradient dephasing. 
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Introduction - The discovery of the very complex vac¬ 
uum in superfluid 3 He proved to be highly stimulating to 
the theory of symmetry breaking [1] and topological de¬ 
fects [2, 3]. New features of quantum gases emerged with 
the first realizations of spinor Bose-Einstein condensates 
(SBEC) [4, 5], thanks to the many degrees of freedom - 
both internal and external - and to the excellent control 
of the experimental parameters. SBECs are extremely 
rich and versatile systems to study complex quantum 
vacuua [6], for example to test the validity of universal 
phenomena like the Kibble-Zurek mechanism [7 9], or to 
study Goldstone modes such as gapless magnons [10]. 

The coupling of SBECs to magnetic fields has been 
exploited to study quantum phase transitions in SBECs 
[11-13], and to realize point-like topological defects such 
as Dirac monopoles [14-16] and 2D skyrmions [17]. In 
these works, spin symmetries and topology were induced 
by strong gradients, e.g. 37 mT/m in Ref. [14]. Here we 
show that via a quantum phase transition at lower gra¬ 
dients ~ 0.5 mT/m, the ferromagnetic interaction can 
“freeze out” the gradient effect. This suggests the possi¬ 
bility of magnetic sensors free from gradient dephasing, 
a practical limitation in coherent magnetometry [18-24]. 

We use group theoretical methods that have proven 
fruitful in the analysis and classification of SBEC phases 
[1, 25]. The interaction with the gradient realizes an in¬ 
teresting spin-position coupling that explicitly breaks the 
parity V and time-reversal T symmetries, only preserv¬ 
ing the combined VT symmetry. Numerically solving 
the Gross-Pitaevskii equations of the system, we observe 
that below a critical value of the magnetic field gradient, 
the VT symmetry is spontaneously broken and a nonzero 
overall magnetization appears. This occurs when the fer¬ 
romagnetic interactions dominate the coupling energy be¬ 
tween the gradient field and the spins, resulting in a glob¬ 
ally polarized condensate. Moreover, we observe that this 
effect is associated with a phase transition. Interestingly, 
discrete, and in particular VT , symmetry breaking is also 
observed in Bose gases with spin-orbit coupling [26-28]. 


System and mean-field energy - We consider a spin-1 
BEC with ferromagnetic interactions in the presence of 
a magnetic field gradient. More specifically, the numer¬ 
ical calculations are performed for a s 'Rb SBEC in the 
F = 1 ground state, which is composed of three Zeeman 
sublevels tof = —, 0, +. 

Within the mean-field approximation, the spin-1 BEC 
is described by a spinorial field with three complex 
components 'P(r) = [^_(r),^> 0 (r), ijj+{ r )] T , where r = 
(. x,y,z) T are the spatial coordinates, and W( r ) * s the 
mean-field wavefunction for the atomic distribution in 
the magnetic sublevel m.F = V- 

The mean-field energy density of the system, coupled 
to a magnetic field distribution B (r), is [29, 30] 
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where the Latin letters designate the spatial coordi¬ 
nates (x, y, z) and the Greek letters the spin coordinates 
(—, 0, +). The first term is the sum of the kinetic energy 
(m is the atom’s mass) and of the trapping potential 
V(r), assumed to be harmonic and spatially isotropic. 
The third term is the energy density resulting from the 
coupling with the magnetic field Bfir), where g is the gy- 
romagnetic ratio, /.is the Bohr magneton, and F l is the 
generator of spin-1 rotations around the i axis. The terms 
containing Co and C 2 < 0 describe the spin-independent 
and ferromagnetic spin-dependent collisional energies, re¬ 
spectively. 

The magnetic field is chosen to be a pure gradient (no 
bias) along the z axis, i.e. the divergenceless field 

B(r) = B\x/ 2, y/ 2, -z) T = B'Ar, (2) 


where we defined the metric A = diag(l/2,1/2, —1). The 
mean-field energy density associated with the gradient 
coupling is thus 


£ g [^]( r ) = ~9TbB' Ar • F(r), 
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where F = (5F X , 3 y , 3 r z ) T and 3/ = V , a-^'a/ 3 V’/ 3 - This inter¬ 
action realizes a spin-position coupling (similar to spin- 
orbit coupling L • F). The total mean held energy related 
to the interaction with the magnetic held gradient is 
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Symmetries - We now study the symmetries of the 
problem, that is, the transformations that leave invari¬ 
ant the mean-held energy £[4*] = f £[4/](r) d 3 r. Fo¬ 
cussing on the invariance of the gradient coupling energy 
[Eq. (4)], we observe that the spin-position coupling ex¬ 
plicitly breaks several symmetries, both continuous and 
discrete. 

We hrst consider the continuous symmetries. In ab¬ 
sence of magnetic held gradient, the energy is invariant 
under both spin and space rotations. When a gradient 
is present, the system is only invariant under combined 
spin-space rotations around the z axis. More precisely, 
we dehne the operator R z (6) acting on 4/ as 

R z (0)^{ r) = e wpz T , (5) 

where L z is the generator of spatial rotations around the 
z axis. The set of all such transformations for 6 € 1R, is a 
one parameter rotation group, denoted as SO(2)p + ^ 2 . 
From Eq. (4), it is straightforward to show that 
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where r' = e~ l6L r. After the change of variable r' —» r, 
we obtain E G [R Z {6 )4>] = Eg[ 4/], proving the invariance 
under In other words, the spin-position 

coupling explicitly breaks the SO(3) F x SO(3) L symme¬ 
try into S0(2) F2+L2 . Therefore, the magnetization F 
and the orbital angular momentum L are no longer in¬ 
dependently conserved, and only the total longitudinal 
angular momentum F z + L z is conserved [31]. 

Beyond continuous symmetries, the mean-field energy 
also exhibits discrete symmetries. We define the spatial 
inversions, corresponding to the parity symmetry, as 


P x : ty(x,y,z) H> (7) 

and similarly for P y and Pr. The spin inversions, corre¬ 
sponding to the time-reversal symmetry, are defined as 
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The energy functional E G is not invariant under parity 
nor time-reversal [32], however it is invariant under a com¬ 
bined space-spin inversion. For example, we consider the 


combined action of P x and T x on the mean-field energy 
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and a change of variable x —» —x leads to £g[P:i-T x 4'] = 
Eg [4']. Therefore, E G is invariant under the discrete 
group VT X = {1,P X T X }. By analogous arguments, the 
mean-field energy is also invariant under VT y and VT Z , 
defined analogously. Therefore, E G is fully symmetric 
under VT = VT X x VT y x VT Z [33]. 

Finally, the energy is also invariant under a phase-shift 
4>(r) —> e l< ^4 , (r), i.e., under the gauge group U(l)^. As 
a consequence, the global symmetry group of the mean- 
field energy is 


G = U(l) 0 x SO(2) F+iI x VT. (12) 

Numerical simulations - Whereas the 11(1)^ symmetry 
is spontaneously broken by the Bose transition, different 
breaking scenarios are possible for the remaining factors 
of the group G, as a result of the competition between 
magnetic field gradient coupling and spin-dependent in¬ 
teractions. Here, using a numerical simulation, we study 
the spontaneously broken symmetries and characterize 
the different phases. 

The mean-field evolution of the system is determined 
by three coupled Gross-Pitaevskii equations associated 
to the mean-field energy density Eq. (1), explicitly 

ihdtip a = [-(fi 2 /2m)V 2 + V(r)] ip a + coip}ippip a 

- g^ B B'Kjr 0 F l a ^ p + c 2 ^*F^ l/ F l a0 ilj fg , 

(13) 

where a € {+, 0, —}. The trapping potential V(r ) is cho¬ 
sen isotropic and harmonic with frequency 100 Hz. The 
atomic species is 8, Rb in the F = 1 hyperfine ground 
state, thus m = 1.44 x 10 -25 kg, g = —1/2, and the colli- 
sional interaction parameters are Co = 5.16 x 10 -51 J-m 3 
and C 2 = —2.39 x 10 - 53 J-m 3 . Since C 2 < 0 the spins 
experience ferromagnetic interactions. 

The stationary state of the GPEs is numerically de¬ 
termined using an imaginary-time method from the 
GPELab 3D solver [34]. The time and space discretiza¬ 
tion is achieved by a backward Euler spectral FFT 
scheme. The initial state is an equal superposition of 
all ?n_p states: 4/(r) = \f f( r) [1,1,1] T , that is, a state 
pointing in the -\-x direction, and where /(r) is a Gaus¬ 
sian distribution with width the harmonic oscillator ra¬ 
dius of the trapping potential V(r). 

As shown on a simulation result in Fig. 1, the magnetic 
field gradient induces a spatial separation of the various 
nriF components along the z axis. We also observe in 
Fig. 1 (b) that the spins in the x — y plane are mainly 
oriented along x, thus the S0(2) F2+i2 symmetry is spon¬ 
taneously broken by the ferromagnetic interaction. This 
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FIG. 1. (Color online) (a) Isoprobability surfaces of the spinor 
components: ip- (yellow, lower), ipo (red, middle) and ip+ 
(blue, upper), (b) Streamlines (curves tangent to the spinor 
field) in three orthogonal planes, and the color map shows 
the total atomic density. The magnetic field gradient is B' = 
3 x 10 -4 T/m and the number of atoms is 5 x 10 5 . 


effect also appears in Fig. 1 (a), where the ip 0 wavefunc- 
tion is not invariant by rotation around z. The prevailing 
direction is x in the current simulations due to the choice 
of initial conditions, however it can be any direction f in 
the x — y plane [35]. 


B' = 0. The sudden change of magnetization when vary¬ 
ing the gradient B' is the signature of a phase transi¬ 
tion, in particular the discontinuity at large atom number 
(N > 10 5 ) would indicate a first order character for the 
transition. However, finite size effects cannot be easily 
excluded from the current simulation and further stud¬ 
ies on the exact nature of the transition should be per¬ 
formed. We call weak gradient (WG) phase the phase 
for B' < B ' c , and strong gradient (SG) phase the one for 
B' > B' c . The phase transition is also visible in Fig. 2: 
R± is independent of A B in the SG phase (A B > 0.7 nT 
for N < 10 5 ). 




FIG. 3. Total magnetization along x versus the gradient B' 

(a) and the magnetic field difference across the condensate 

(b) . Same legend as in Fig. 2. 
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A B = R±B' (nT) 


FIG. 2. Cloud radius over the Thomas-Fermi radius versus 
B for various atom numbers. The legend at the top indicates 
the number of atoms. 

A natural measure of gradient strength is A B = R±B ', 
where R± is the mean-square radius of the ± compo¬ 
nents. As shown in Fig. 2, this radius, normalized by the 
Thomas-Fermi radius Rtf, is nearly independent of the 
atom number in the range 2000 < TV < 10 s . Below this 
range the low density voids the Thomas-Fermi approxi¬ 
mation, and above this range the large size of the system 
exhibits the first-order behaviour of the transition (see 
next paragraph). 

In Fig. 3 (a), we show the behaviour of the total mag¬ 
netization along x, (F x ) = f J x (r)d 3 r [36]. We observe 
that, below a critical value of the gradient B' c of the or¬ 
der of 0.5 mT/m, the cloud is spontaneously magnetized 
along x, converging towards a fully polarized state for 


Broken symmetries - From the simulation results, we 
can study the WG and SG phases from a symme¬ 
try perspective. As already pointed out, U(l)^ and 
SO(2) Fii+LI are spontaneously broken as a result of 
the Bose/ferromagnetic transition. Thus, all continu¬ 
ous symmetries of G are broken, and only discrete sym¬ 
metries remain. In Fig. 4, we compare the spinor field 
streamlines, within slices of three orthogonal planes, for 
condensates in both phases. As pointed out earlier, here 
the prevailing axis f is oriented along x. 

In the WG phase [Figs. 4 (a) and (b)], the system is 
not invariant under P f Tj, that is, the space-spin inversion 
along the axis f. Conversely in the SG phase [Fig. 4 (c)], 
the system is invariant under both Pj-T^ and Pj Tj , 
where fj_ is orthogonal to f and lies in the x — y plane. 

The symmetries of the order parameter in a given 
phase constitute the isotropy group , which is a sub¬ 
group of G. For the WG phase, the isotropy group is 
%g = VTf x VT Z , where VT} = {l, PjTj}, whereas 
for the SG phase it is Hsg = VT $ x VTf ± x VT Z - 
Order parameter space - Many properties of a given 
phase can be understood from the broken symmetries 
[2], mathematically defined as the quotient of the overall 
symmetry group by the isotropy group H of the phase: 
R = G/H, called the order parameter space. If is 
an arbitrary order parameter of a given phase, called the 
standard order parameter , then the order parameter man¬ 
ifold A4 of the phase results from the action of R onto 
\k 0 , that is, A4 = Rty 0 [37]. It results that many crucial 
properties of a phase, fully defined by the manifold A4, 
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FIG. 4. (Color online) Spin streamlines (solid lines with ar¬ 
row) and total density (color map) for N = 2 x 10°. Slices 
in three orthogonal planes are presented, from left to right 
are slices in the y — z, x — z and x — y planes, respectively, 
(a) WG phase with B' = 0.4 mT/m, (b) WG phase with 
B' = 0.535 mT/m, and (c) SG phase with B' = 0.6 mT/m. 


are embedded in the order parameter space R. This is in 
particular true for the topological properties. 

The order parameter space of the WG phase is [38] 

Rwg = G/Hwg = U(l)^ x S0(2) F2+L2 . (14) 

From a topological perspective, both U(l)^, and 
S0(2) F2+L2 have the topology of a circle ffi 1 , and the 
order parameter space is a torus: i?wG - S 1 x S 1 = T 2 . 
For the SG phase, the order parameter space is 

Rsg = U(l) 0 x (S0(2) F2+i2 /VTt) • (15) 

Let 0~ f be the angle between the vector f and the x axis, 
then the quotient by the group VT$ identifies the ele¬ 
ments of S0(2) F2 , L z according to the equivalence rela¬ 
tion ~, defined as for all 0, O' £ R, 

R z (9) ~ R z {0’) 4=> (O' = 0 or O’ = 20 j. — 0 ) mod 27r. 

(16) 

As a consequence, the circle associated to S0(2) F2+i2 
becomes a closed line, topologically equivalent to a closed 
interval: SO(2) pz+LZ /VTf — I, and the order parameter 
space of the SG phase is a cylinder: Rsg — S 1 x I. 


The knowledge of the order parameter space topol¬ 
ogy provides information on topological defects, in par¬ 
ticular their stability. For both phases, the order pa¬ 
rameter space is connected and therefore domain walls 
are unstable. However, these spaces are not simply 
connected, and thus stable ID defects such as vortices 
can form. In the SG phase the fundamental group of 
Rsg is 7Ti (S 1 x I) = Z, and vortices are classified by 
a single winding number. Whereas for the WG phase, 
7Ti (T 2 ) = Z x Z and the vortices are characterized by 
pairs of winding numbers, one related to the superfluid 
phase and the other to the magnetization. Moreover, 
higher order homotopy groups are trivial, and in partic¬ 
ular point-like topological defects are unstable. 
Conclusion and outlook - We have studied the ground- 
state properties of a ferromagnetic F = 1 spinor conden¬ 
sate in the presence of a magnetic field gradient, a con¬ 
figuration that breaks both time-reversal T and parity 
V symmetries, but preserves the combined VT symme¬ 
try. Simulation reveals a phase transition that sponta¬ 
neously breaks also this VT symmetry for weak gradi¬ 
ent strength. Distinct topological defects are predicted 
in the weak- and strong-gradient phases. The fact that 
the polarization of the WG phase, parallel to f, is free 
to precess about the z axis while protected by a phase 
transition suggests an attractive system for coherent field 
sensing. In contrast to other atomic field sensors [18-24], 
gradient-induced dephasing may be frozen out by the fer¬ 
romagnetic interaction. This possibility motivates the 
study of the dynamics of this system under a combined 
gradient and bias field. 
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